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RESUMEN 

Presentamos resultados de los calculos para HD 5980 efectuados con nuestro modelo de interaccion por fuerzas 
de marea. El modelo predice una mayor actividad atmosferica despues del paso por el periastro, asi como una 
mayor actividad ocasional en zonas intermedias entre el ecuador y el polo. 

ABSTRACT 

The LBV/WR system HD 5980 contains a short-period, eccentric binary system with interacting stellar winds. 
In this paper we summarize results from model calculations of the tidal flows on the LBV component showing 
that energy dissipation rates, E, associated with turbulent viscosity are orbital-phase dependent as well as 
variable over the stellar surface. We speculate that if E contributes towards driving mass-loss, the strongest 
wind-wind interaction effects may occur after periastron passage. In addition, the model suggests the presence 
of stronger outflows localized at polar angles -50° during part of the orbital cycle. Thus, the analysis 
of wind-wind interactions in this system requires that models be revised to incorporate non-stationary and 
asymmetric wind structures. 
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1. ERUPTIONS, WIND- WIND INTERACTIONS 
AND TIDES 

HD 5980 is an amazing system. Niemcla (1988) 
was the first to point out that it consists of two Wolf- 
Rayet-like components referred to as star A and star 
5 in a relatively close, eclipsing and excentric orbit 
(P=19.3, e~0.3) and a third source, referred to as 
star C, that may simply lie along the line-of-sight to 
the close pair. The spectral characteristics and vi- 
sual brightness of the system underwent significant 
changes between the late 1970's and 1993, when it 
entered an eruptive state that lasted ~1 year (Bate- 
son & Jones 1994; Barba et al. 1995; Koenigsberger 
et al. 1995). The activity involved an increase in vi- 
sual brightness and mass-loss rate, and a decrease in 
wind velocity and effective temperature, similar to 
the eruption phenomena observed in luminous blue 
variables (LBVs). The radial velocity variations ob- 
served in the very rich emission-line spectrum that 
appeared after the eruption led Barba et al. (1996; 
1997) to conclude that the instability producing the 
outbursts originated in star A. A detailed review of 
HD 5980's properties is provided by Koenigsberger 
(2004). 
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The ZAMS masses of star A and star B are in- 
ferred to be >100 Mq (Koenigsberger 2004). The 
substantial mass loss required for them to have 
reached their current masses (M^ ~50 Mq, Mb ~30 
Mq, Niemela et al. 1997) may have been achieved 
through multiple events as those of 1993/1994. 
LBVs are associated with the evolutionary state 
during which large quantities of mass are ejected al- 
lowing the star to reach the W-R phases with highly 
depleted hydrogen envelopes. With the possible ex- 
ception of 77 Carinae, there is no known LBV with 
such a developed W-R spectrum as that displayed 
by HD 5980. 

1.1. A changing wind-momentum ratio 

The spectrum of HD 5980 in the late 1970's, 
with its broad He II and N V lines, was typical 
of the "early" W-R stars of the nitrogen sequence 
(WNE; van dcr Hucht, 2001). This spectrum is be- 
lieved to originate in the wind of star B. Over the 
next decade, however, numerous lines from lower- 
ionization atomic species appeared, implying a grow- 
ing presence of a cooler stellar wind, which we now 
attribute to star A. Clearly, the emerging dominance 
of star A's wind implies changes in the wind- wind in- 
teraction (WWI) region characteristics. 

Emission-line profile variations observed in opti- 
cal and UV wavebands are phase-locked and should, 
in principle, provide information on the geometry of 
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velocity (-150 km s"*) 

Fig. 1. Line profiles of tlie semi- forbidden N IV emis- 
sion line at A 1486 A observed at two different orbital 
phases (0 ~0.40 and 0.83 -dots) for each of three differ- 
ent epochs (1991, 1995 and 1999). The qualitative na- 
ture of the line-profile variations remains constant. Sets 
of profiles for different epochs are vertically displaced for 
clarity in the figure. 

the changing WWI region (Moffat et al. 1998). Sur- 
prisingly, however, the nature of the variability has 
remained the same ever since it was first reported 
by Breysacher & Westerlund (1978) and quantified 
by Breysacher, Moffat & Niemela (1982). The varia- 
tions consist of periodic changes in width and degree 
of asymmetry as a function of orbital phase. Figure 
1 illustrates the NIV] 1486 A emission line variabil- 
ity: it is always narrower and sharply peaked near 
the eclipse when star B is "in front" ((/) ~0.40), while 
becoming broader and weaker when both stars are 
unocculted at (j) =0.83. The three epochs that are 
displayed correspond to prc-eruption (1991), post- 
eruption (1999) and ^ 1 year after maximum (1995). 
This persistent trend is unexpected since the geom- 
etry of the WWI regions depends on the momen- 
tum ratio of the stellar winds, a ratio that changed 
over time as star A's wind became more dominant. 
Thus, HD 5980 appears to provide another example 
of the discrepancies that arise when confronting the 
current WWI models with observations (sec presen- 
tations in this Volume by Rauw, Pollock, Williams, 
among others). But at the same time, because of the 
large observational database available for HD 5980, 
its behavior may provide a clue to identifying the 
source of the discrepancies. 



1.2. Tides and non- stationary, asymmetric winds 

Like many of the intriguingly active binaries, 
HD 5980 has an excentric orbit. In such systems, 
the tidal forces are time- variable. The preliminary 
exploration of tidal effects in HD 5980 led to the 
conclusion that they are non-negligible (Koenigs- 
berger et al. 2002), and thus raised the question 
of whether they may be responsible for some of the 
system's peculiarities. Recent results of our calcula- 
tions (Moreno & Koenigsberger 2007) indicate that 
the tidal fiows near the stellar surface can liberate 
considerable amounts of energy through dissipative 
processes. The magnitude of the energy dissipation 
rate, E, depends on the stellar and orbital param- 
eters. In the following sections, we'll describe the 
two basic conclusions of the calculations: 1) maxi- 
mum E occurs after periastron passage, not at peri- 
astron; and 2) at certain orbital phases, larger val- 
ues of E are generated at intermediate polar angles 
than in the equatorial belt. These results are rel- 
evant for WWI theory if we assume that E con- 
tributes towards enhancing the stellar mass-loss rate. 
The winds in excentric binaries such as HD 5980 
would then be intrinsically non-spherically symmet- 
ric and time-dependent, thus leading to discrepan- 
cies when comparing the observational diagnostics 
of WWI with the predictions of stationary models. 

2. ENERGY DISSIPATION FROM TIDAL 
FLOWS IN HD 5980 

Tidal effects are important when a star's rota- 
tion, Lu, is not synchronized with the orbital motion, 
57. In excentric binaries, this is generally the case 
since the orbital motion varies with phase while the 
stellar rotation rate remains constant. There is no 
direct observation of vsini for star A or star B. But 
the low-amplitude variations of star C's narrow ab- 
sorption lines over the 19.3-day orbit has been in- 
terpreted in terms of contamination by very broad 
absorptions arising in star A but not visible due to 
their superposition on the emission lines (Georgiev 
& Koenigsberger 2004). Under this assumption, 
vsini >200 km/s, thus providing an estimate for the 
ratio uj/flper ^^2.33, where Qper is the orbital angu- 
lar velocity at periastron. This means that star A 
rotates super-synchronously throughout the orbital 
cycle. 

The basic method is described in Moreno & 
Koenigsberger (1999) and Moreno et al. (2003). 
With the recent extensions (Moreno & Koenigs- 
berger, 2007), the code now computes the ampli- 
tudes of the tidal flow in a thin surface layer over the 
entire stellar surface. These amplitudes are used to 
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Fig. 2. Predicted energy dissipation rate per unit density 
due to tire tidal flows on the stellar surface of star A as 
a function of orbital phase. Maximum rates occur after 
periastron passage. Star A is "in front" at <j) =Q, and 
periastron passage occurs at (f) ~0.07. 

estimate the shear energy dissipation rates, E that 
arise from the relative motions of different surface 
layers using an extension of the approach described 
in Toledano et al. (2006). 

Figure 2 illustrates the time-dependence of E as 
a function of orbital phase from our HD 5980 model 
calculation. The first thing to note is that the max- 
imum rates are generated after periastron passage, 
(j) ~0.l-0.25, with minimum values around (f) =0.8. 
Thus, it is tempting to suggest that the persistent 
line-profile variability shown in Figure 1 may be as- 
sociated with the orbital-phase dependent changes 
in E. But in what way do the tidal effects produce 
these variations ? 

The answer may lie in the distribution in latitude 
of E. Figure 3 illustrates E = E{9) for several dif- 
ferent times within the orbital cycle. At periastron 
(dotted curve) and until ~2 days thereafter, maxi- 
mum E{9) occurs at the equator and systematically 
decreases towards the pole {9 =0 at the pole). By 
day ^^3 and until apastron, however, there is a dis- 
tinct change in this trend whereby maximum E{9) 
now occurs at 9 ~30-50°. If wc assume that E ^ M, 
the results shown in Figure 3 imply that between 
days ~3 and ~6 after periastron passage (0 ~0.2 - 
0.4) mass outflow from regions at intermediate po- 
lar angles is more intense than from the equator. A 
denser "polar" wind would produce a narrower and 




Latitude (degrees from pole) 

Fig. 3. Predicted energy dissipation rate per unit density 
due to the tidal flows on the stellar surface of star A as 
a function of latitude at periastron and several times (in 
days) after periastron passage. Note that E is greater 
at intermetidate polar angles {6 ~20-50 °) than at the 
equator around ~days 3-6 after periastron. 

sharply peaked emission-line profile, as seen in Fig- 
ure 1 for (p ~0.40. Additional observational evidence 
for polar outflows within this orbital phase interval 
is provided by Villar-Sbafh ct al. (2003) who state 
that "the mass-loss of HD 5980 around cj) =0.36 pre- 
sented fluctuations in axial symmetry ranging from 
very rapid density enhancements along the orbital 
plane to polar ejections." At orbital phase cj) ^0.8, 
E{9) has a relatively small gradient between the pole 
and the equator and is significantly weaker than near 
periastron. Thus, the stellar wind structure of star A 
should be more spherically symmetric at this phase. 

Clearly, there are additional contributions to pro- 
ducing the line-profile variations, such as the physi- 
cal and wind eclipses, as well as contributions from 
the WWI zone and all of these need to be consid- 
ered. However, it is encouraging that the asymme- 
tries in wind structure predicted by the tidal inter- 
action model arc consistent with the emission-line 
profile variability. 

3. FINAL REFLECTIONS 

A one-layer tidal interaction model for HD 5980 
predicts that the energy dissipation rate due to the 
tidal shearing flows are time-dependent and non- 
spherically symmetric. Speculating that E con- 
tributes towards the stellar wind mass-loss, leads 
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to the conclusion that M may also be locally en- 
hanced near specific surface locations. In particu- 
lar, outflows at intermediate polar angles may be 
stronger than at the equator at particular orbital 
phases. Furthermore, the model predicts that max- 
imum E should occur after periastron passage, thus 
implying an overall enhanced M compared to other 
orbital phases. 

Post-periastron events seem to occur in a wide 
variety of binary systems, such as WR 140, jy Cari- 
nae and others, raising the question of whether the 
tidal effects described above are more prevalent than 
one may have anticipated. Is it possible that the 
stronger WWI effects ("outbursts") that occur af- 
ter periastron passage are associated with stronger 
mass-loss rates at these phases induced by tidal in- 
stabilities ? If this were the case, the source of the 
discrepancies between wind-wind interaction model 
predictions and the observations may simply reside 
in the assumption of stationary and spherically sym- 
metric winds. 

It is interesting to note that the hypothesis of 
enhanced M arising from E may not be entirely un- 
reasonable. Given HD 5980's huge UV luminosity 
(see, for example, Koenigsberger 2004), it is likely 
to be on the verge of the Eddington Limit, as other 
W-R stars appear to be (Goeffncr, this workshop). 
Thus, small additions of energy in sub-photospheric 
layers could drive it to a super-Eddington state. We 
speculate that viscous shear energy dissipation re- 
sulting from the tidal forces may be a non-negligible 
contributor to this small needed additional energy. 

Within this context, a final consideration con- 
cerns the sudden eruptive events in HD 5980. Mon- 
itoring of HD 5980 at visual and UV wavebands 
prior to, during and after its eruptions has yielded 
a unique data set that provide clues for constraining 
the eruption mechanism. For example, ultraviolet 
observations suggest that the onset of the eruptive 
state involved rapid transitions between a fast and a 
slow stellar wind (Koenigsberger 2004). Hence, it is 
likely that the eruption occurred when the wind be- 
came so dense that the bistability limit was crossed 
(Lamers, Snow & Lindholm 1995). But this is only 
the symptom of a more deep-seated phenomenon 
that causes the instability leading to the enhanced 
density wind in the first place. Tidal effects are very 
sensitive to the star's radius. If HD 5980 (and other 
similar stars) are undergoing an evolutionary tran- 
sition by which outer layers are expanding, the am- 
plitudes of the tidal flows are expected to grow sig- 
nificantly. If the hypothesis that E — > M can be 
shown to stand on firm ground, this would provide 



a mechanism to remove significant amounts of mass 
as the star tries to evolve towards the red end of 
the H-R Diagram. Whether the mass-shedding oc- 
curs as episodic eruptions, such as we've observed 
in HD 5980, or whether it is through a sustained 
high-density wind, requires an understanding of the 
E ^ M process. Since our model neglects the effects 
of intrinsic stellar oscillation modes, effective tem- 
perature variations and radiation pressure, we are 
unable to go beyond the speculative realm at this 
time. 

GK thanks Jeff Kuhn and Stan Owocki for very 
helpful discussions. Support from PAPIIT/DGAPA 
grant INI 19205 is gratefully acknowledged. 
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Transcription of questions and answers 
Moffat: There are other sources of perturbation 
that should be considered, such as the heating effects 
from the X-rays produced at the bow-shock head be- 
tween stars A and B, as we started to do in Moffat 
et al. (1998). 

Reply: Wind-wind collisions and the effects associ- 
ated with them are inevitable if the winds of both 
stars are able to accelerate to large speeds in the 
intervening region. However, strong tidal perturba- 
tions are also inevitable if the stellar rotation is not 
synchronized with the orbital angular velocity. 
Nathan Smith: The high luminosity high mass loss 
make LBV ourbrust and the fact that its in a binary 
naturally make me think of comparing HD 5980 to 
Eta Car. In that case, its curious that their outbursts 
were so different. Namely, Eta Car ejected over 10 
Mo in its 19th century outburst, whereas HD 5980 
only ejected about 0.001 Mo ...that's four orders of 
magnitude different and it begs for an explanation, 
given that the primary star's luminosities and mass 
loss rates are within a factor of 2. One obvious dif- 
ference that comes to mind is that Eta Car has a 
5 year period whereas HD 5980 has amuch smaller 
separation with only a 19.3 day period. So I am won- 
dering if the closer companion in HD 5980 somehow 
regulates the LBV instability so that it erupts more 
often with less mass, and therefore doesn't build up 
the catastrophic instability that leads to a major out- 
burst like Eta Car. Any thoughts on that ? 
Reply: We do not know what the underlying mech- 
anism for the instability in HD 5980 really is. If we 
assume that its radius is growing due to evolution- 
ary processes, it is possible that the presence of the 
binary companion produces a "premature" eruption 
through the tidal oscillations mechanism, removing 
mass, and thus slowing (momentarily) the expan- 
sion. The closer the binary companion is, the more 
frequent such events would be expected to occur. 
MoflFat: (After N. Smith comparing HD5980 to 
Eta): Eta at periastron could be more similar to 
HD 5980, given its large eccentricity. 
Reply: The orbital separation in HD 5980 is only 
-100 Ro. 

Maeder: What about the chemical abundances 
which may tell us something on the evolutionary 
stage of HD 5980 ? 

Reply: The Non-LTE analysis that we made on 
the eruptor's spectrum (Koenigsberger et al. 1998) 
indicates that there is a significant fraction of hy- 
drogen. A determination of other chemical element 
abundances would be highly desirable. 
Walborn: Could the star be hitting the Eddington 



Limit as well as or instead of the Bistability Limit 
as it tries to evolve ? If there are internal magnetic 
fields combined with differential rotation and tidal 
oscillations, the effects may also be catastrophic. 
Reply: Yes, indeed. 

Stan Owocki: (replying to Nolan): One idea for a 
way to trigger a super-Eddington luminosity is that 
the tidal excitation of pulsations you mentioned is 
associated with a breaking of the spherical symme- 
try that is essential for blocking the radiation in the 
envelope. The associated clumping or "porosity" of 
the envelope could then allow radiation to escape 
from an edge, leading to super-Eddington brighten- 
ing that could then drive the mass eruption. 
Nathan Smith: (after Stan's comment): Did the 
bolometric luminosity of HD 5980 change during its 
LBV outburst ? Pinning that down is important 
for explaining the differences between this event and 
Eta Car. If HD 5980 stays at roughly constant L^oi , 
the likely implication is that it is crossing the bista- 
bility jump and develops a pseudophotosphere as in 
a normal S Dor outburst, whereas Eta Car violated 
the classical Eddington limit and therefore suffered 
a much more violent, deep seated event. In that 
context, one question is does HD 5980 show any ev- 
idence for a massive circumstellar shell ejected in an 
ancient EtaCar-like eruption ? 

Reply: We (Koenigsberger et al. 1998, ApJ 499, 
889) derived Lboi =3x10^ Lo, but Drissen et al. 
(2001, ApJ 545, 484) derived --10^, so it is not clear 
whether there was a change or not. 
Moffat: Tidal oscillations may indeed be an im- 
portant trigger of the eruption of star A. But once 
star A reaches maximum size, it is much larger than 
the orbital separation. This can lead to common en- 
velope evolution for a time, and eventually shorten 
the orbital period leading in the (very) short period 
WR4-0 binaries we see in the Magellanic Clouds. 
This means that the LBV phase may indeed be cru- 
cial for explaining WR stars, given the reduced mass- 
loss rates of their progenitor 0-stars. 
Reply: The common envelope phase (if we can call 
it that!) during the 1994 eruption lasted less than a 
year, probably too short to cause any changes in the 
orbit. But if this were to happen frequently enough... 



